QCD Saturation Equations including Dipole-Dipole Correlations 
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We derive two coupled non-linear evolution equations corresponding to the truncation of the 
Balitsky infinite hierarchy of saturation equations after inclusion of dipole-dipole correlations, i.e. 
one step beyond the Balitsky-Kovchegov (BK) equation. We exhibit an exact solution for maximal 
correlation which still satisfies the same asymptotic geometric scaling as BK but with the S-matrix 
going to 1/2 (instead of 0) in the full saturation region. 

1. In perturbative QCD, parton saturation, i.e. the modification of the distribution of quarks and gluon distributions 
in a target, is known to lead to an infinite set of coupled evolution equations in energy for the correlation functions of 
multiple Wilson lines 0. This set of equations is expected to be equivalent to the Jalilian-Marian, Iancu, McLerran, 
Weigert, Leonidov and Kovner (JIMWLK) functional equation 0- In the approximation where correlation functions 
for more than two Wilson lines factorize, i.e. (tr U 2 tr U 2 ) ~ (tr U 2 ) (tr U 2 ), the problem reduces to an unique non- 
linear Balitsky-Kovchegov (BK) equation jlj, 3] for the dipole density. This amounts to considering a large- N c limit of 
independent dipole-target collisions, e.g. on a large nucleus. Recently it was shown that the translation-invariant (no 
impact parameter 6-dependence) non-linear BK equation lies in the universality class of the Fisher and Kolmogorov- 
Petrovsky-Piscounov (F-KPP) equation |^|, leading to asymptotic traveling wave solutions in the transition region to 
saturation. This provides a mathematical realization of the phenomenologically motivated j5[ geometrical scaling. In 
addition, the BK equation leads to an S-matrix element which goes to in the full saturation region. Note that, if 
one assumes that translation invariance is a good approximation for scattering near 6 = 0, this amounts to the black 
disk limit a e i/a tot {b ~ 0) — >• 1/2. 

The problem we want to address is the study of QCD saturation including the effects of nontrivial two-dipole 
correlations. To this end we have to proceed one step further in the hierarchy of Balitsky equations [1|. Our main 
observation is that one can perform a truncation of this set at the level of two-dipole correlations, as follows: 

• We keep only dipole-like terms (i.e. we neglect higher multipoint traces like (tr C/ 4 ) and (tr U 6 )). 

• We keep full two-dipole correlations and only neglect independent three- and higher dipole correlations. 

In this paper we derive the general closed set of two non-linear equations for the dipole density and the dipole- 
dipole correlation function, see equations H9I10JI below. Moreover, using the translational invariance hypothesis (no 
6-dependence) , we find a particularly convenient form of these equations in a Fourier transform representation see 
equations (|13I14|) . We then exhibit a solution in terms of a modified BK equation when assuming maximal correlation, 

1. e. when the correlation stays independent of the separation distance between the two dipoles see <|21I22|) . We analyze 
the solutions, examine the modifications with respect to BK ones and discuss the prospects of our study. As an aside 
we note that some hierarchies of evolution equations have been studied in a statistical physics context in 

2. Let us briefly recall the standard derivation of the BK equation in the Balitsky's framework 0. One defines the 
dipole operator 

n, ir r,.r; . (l) 

where Xij are the transverse coordinates of the end point quark and antiquark of a QCD dipole. One calculates its 
evolution with rapidity d/dY using contractions between pairs of U's and virtual corrections for each U (see e.g. eqs. 
(119,120) in 0].) 
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Within the infinite set of saturation equations, the one for the evolution of the dipole operator Q expectation value 
with rapidity reads: 



^7(^01} = £3 / </-.<-2 ca.2-p->i) - -v, .Pn,;; fv-„,j (2) 



where the kernel -K021 is 



^021 — —3 — 2~ ' (3) 
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i.e. the real dipole splitting part in the BFKL kernel |g. Let us now decompose the expectation values into 1-point 
functions and connected correlation functions through 

(D02) = d 02 (4) 
(V Q2 V 21 ) - (V Q2 ) (V 2 i) = d 0221 (5) 

The BK equation is obtained by neglecting do 22 \, i.e. the connected part of the two-dipole correlation function, with 
two coincident points at x 2 . 

Our first goal is to keep now this term, and derive an equation for this quantity. In order to obtain a well-defined 
closed set of equations, it turns out to be necessary to write an equation for the correlation of arbitrary two dipoles, 
not constrained to have one point in common. We thus consider more generally d {T>q 2 T> 2 >i) /dY and perform all 
contractions between the unitary matrices. 

For the truncation, we proceed as follows. All single (1-point) and double contractions within the same dipole 
give essentially two copies of © with the other dipole included as a spectator. To this we have to add contractions 
between U's belonging to different dipoles. These generate only non-dipole terms of the form (tr U 6 ) and (tr U 4 ) 
which we neglect by assumption. The resulting equation reads 



_d_ 

dY 



(DatVvi) = ^3 / d2 ^3 [(^o 3 ^32^2'i) - N c (VmVvi)] K m2 + [{V 2 ,zD sl V Q2 ) - N c (P 02 £>2'i}] K 13V . (6) 



When 2 = 2' this reduces to the expression in [l| up to the neglected (tr t/ 6 ) and (tr C/ 4 ) terms. Keeping track of all 
two-dipole correlations, the expectation values are evaluated giving 

\D03D 32D 2'l) — ^03^32^2'! + ^0332^2'! + e?322'lrf()3 + rfo32'1^32 ■ (7) 

It is important to notice the appearance here of two-dipole correlation functions with four non coinciding points even 
if we would start with 2 = 2', due to the contractions between the first and third T> operators. This justifies the need 
to get the corresponding 4-point evolution equations. Indeed, neglecting this term would lead to spurious divergences 
in the equation. In this way we get: 



^ f 

7^7^022-1 = g-^ / d 2 x 3 {d 03 d 322/1 + d 32 d 03 2'i - N c d 022 > i} K 032 + {(0, 1) (2', 2)} 



(8) 



8 N oi = ^T~T I d ' 2x 2 i N 02 + N 2l - N 01 - N 02 N 21 - N 022 i] K 02 i . ( 1.0 ) 



Going from the Wilson line correlation traces d to dipole densities (dot = N c (l — Nqi) and do 2 2'i = N% ^Vo22'i) w e 
get the final closed set of equations: 

d q 2 N r 

^^022'! = V^T / ePx3{[N322'l+No32'l-N 22'l]K 032 - [N 3N 32 2'l+N 32 No32'l] K 032 } + {(0,1) ^ (2',2)}(9) 

3. It is well known 3] that after imposing translational invariance (i.e. restricting to impact-parameter independent 
regime) and going over to momentum space the BK equation greatly simplifies. We will perform an analogous 
procedure for our set of equations (|§|)- (|10fl . Let us introduce the following 2-dimensional Fourier transfoms 

N 01 = x 2 01 j fh.e ik ^ N k (11) 

AWi = Zq 2 4i / ^^d 2 Q e 1 ^ e I9 '-^'i e iQ < x ™' +x ^ N qq , >Q (12) 

/ Z7T Z7T 



Here q and q' are 2-vectors conjugate to the dipole vectors, while Q is the variable dual to the 2-vector drawn between 
the center of masses of the two dipoles. The nonlinear terms now become local in momentum space, while the linear 
terms reduce essentially to the BFKL kernel [9j. We get finally: 

■Hy N * = Vx(~d L )N q -N 2 q -J d 2 Q JVq )( ^),q} (13) 

d q 2 N 

QY N **,Q = l^T { 2 x(-d L )N qql>Q + 2 X {-d L ,)N qq , <Q - (N q+Q + N q _ Q + N q , +Q + N q -- Q )N qq ,, Q } , (14) 
where L = log(q), 

X(l) = 2^(1) - ^(7) - 1>(1 ~ 7) (15) 

and x (—9l) is an integro-differential operator which may be defined 4] with the help of a formal series expansion 
around some given 70 between and 1, i.e. for the principal branch of the function x- Note that the integral term in 
the first equation (fT3)l corresponds to the integration of a dipole-dipole correlator (5) with the kernel K 2i over the 
separation distance of the two dipoles with one coinciding endpoint. 

The initial conditions for correlations can be formulated in the variable Q by assuming a gaussian form e <=<"-»-. ^ 
corresponding to a transverse correlation length l cor r.- In order to illustrate the impact of these correlations on the 
problem, let us consider the limiting case of maximal correlations i.e. when the dipole-dipole correlations are local 

in Q (lcorr. —* Oo)'- 

N qq >, Q =N qq ,5 2 (Q) . (16) 

Then the equations simplify and we get 

w Nq = ^~ 9l ^ - N " - ^ 

d a 2 N 

OY N m> = Vx(-d L )N qq/ + 2 X {-d v )N qql - 2{N q + N q ,)N qq ,} . (18) 

Remarkably enough this set of equations admits an exact solution (apart from the trivial solution N qq : = 0) using 
the ansatz N qq i = \N q N q >. Note that this form may be quite plausible physically as typically these dipoles are 
separated by very large distances, therefore we do not expect strong dependence of the correlations on the relative 
sizes. Plugging the ansatz into equations Ijl7|l - H18|l we obtain A = 1 and find the following system 

N qq > = N q N q , (19) 
F) n 2 N 

W N " = { 2 *HW - 2iV g 2 } , (20) 



or equivalently in position space: 



iWl = ^02^2-1 (21) 

d a 2 N [ 

ffy N oi = / d2x2 t^ 02 + ^1 - N oi ~ 2N 02N 2 i] K 021 . (22) 

In fact this solution can be extended to the full 6-dependent set of equations (|9I1U|I . 

It is important to note the factor 2 in front of the non-linear term which is the consequence of the correlations. Note 
that a modification of the BK equation of the type (|22|l but with 2 replaced by some constant has been proposed in 
|l2l | on phenomenological grounds to account for the effect of correlations in nuclei. Our approach gives a derivation 
of such a modification from the JIMWLK-Balitsky framework in the special limit of infinite range correlations and 
opens up a way to study corrections due to finite correlation length. 

The solution of the equation l|22() can be obtained from the solution of the uncorrelated BK equation through the 
relation 

N Q1 (l corr . ^00) = ^Nrf x K . (23) 

The S matrix is consequently 5=1 — N = 1 — ^N BK and satisfies the equation: 

3 a 2 N f 

g^Soi = ^Tj d2x * [2(^2 - I) x (Sai - I) - (S01 - *)] ^021 . (24) 



A consequence of the above equations is that the solution for the case with maximal correlations leads to a saturation 
regime where the S-matrix goes to \ instead of 0. Physically, it means that the ratio of the elastic over total cross- 
section (at small impact parameter at least) a e i/atot(b ~ 0) — > \ instead of j for the BK equation. It is remarkable 
that the fixed point solution S — ► J saturates the so-called Pumplin bound 11], which states that a e i + add < \&tot-, 
where add is the contribution of inelastic diffractive channels to the total cross-section. Indeed, a simple calculation 
gives a e i = add =\otot = \- It is interesting to notice that the phenomenological extraction of dipole-proton S-matrix 
values seems to be consistent with S >\. 

Another interesting aspect is that geometrical scaling and more generally, the transition to saturation is not expected 
to be modified, since it follows directly from the exact relation Ij23|l . It is yet another consequence of the universality 
properties of the traveling wave solutions Q. The expectations from the mathematical properties of non-linear 
equations discussed in 4] show that the behaviour of the solution in that region are universal, i.e. independent of 
the precise form of the non-linear damping. Indeed, since equations H13fl - I|14[) have a similar linear term as in the BK 
equation and lead to geometrical scaling both in the maximal correlation regime and in the regime of no correlations 
at all, therefore it is tempting to conjecture that geometrical scaling holds even when the assumption of maximal 
correlations is relaxed. This point certainly deserves further study. 

4. In conclusion, we have shown that a consistent truncation can be performed on the infinite hierarchy []J of Balitsky 
equations for n-point correlators of Wilson lines by keeping track of nontrivial two-dipole correlations. It results in 
a system of two coupled non-linear equations. It is interesting to note that the assumption of maximal correlations 
does not modify the traveling wave picture of geometrical scaling. In contrast, the fully saturated region is strongly 
modified, since the S-matrix admits a different limit. 

Recently, it has been argued that a certain treatment of fluctuations beyond the BK equation may lead to deep 
modifications including strong geometrical scaling violations (l4|. We do not find this phenomenon when including 
dipolc-dipole correlators with maximal (very long-range) correlations. It would be interesting to look what is the 
behaviour of the solution for finite range correlations. A numerical simulation of our system of equations (| 1 3 1> — (| 1 4|) 
seems feasible and could be helpful. 

An important point to study is the dependence of the solutions on the initial conditions, in particular the stability 
of the fixed point solution for maximal correlations when one starts from generic initial conditions (apart from the 
initial condition with no correlations which trivially reduces to the standard BK equation). Interestingly enough it 
seems not out of reach to generalize the maximal correlation ansatz (|21f) to multipoint correlations in the JIMWLK 
context. 

It is interesting to ask the question whether the connected dipole correlation functions could be directly related 
to physical observables. For instance, Fourier transforming back the function N qq i : Q (see II14J) over q,q' leads to a 
relation with the interaction amplitude of 2 dipoles with the target. By crossing symmetry 1 1 "J . some information can 
be obtained on the semi-inclusive dipole(r) — > dipole(r') + X scattering with transverse momentum Q on the target. 
This problem deserves more study, since it could give an interesting link with inelastic diffractive processes. 
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